Lattice two-body problem with arbitrary finite range interactions. 
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We study the exact solution of the two-body problem on a tight-binding one-dimensional lattice, 
with pairwise interaction potentials which have an arbitrary but finite range. We show how to 
obtain the full spectrum, the bound and scattering states and the "low-energy" solutions by very 
efficient and easy-to-implement numerical means. All bound states are proven to be characterized 
by roots of a polynomial whose degree depends linearly on the range of the potential, and we 
discuss the connections between the number of bound states and the scattering lengths. "Low- 
energy" resonances can be located with great precission with the methods we introduce. Further 
generalizations to include more exotic interactions are also discussed. 
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I. INTRODUCTION 

The study of quantum mechanics in periodic structures 
is one of the central topics in condensed matter physics 
since many decades The behavior of electrons in a 
crystal or, more generally, interacting particles in a peri- 
odic potential, even at the few-body level, puts forward a 
major theoretical and numerical challenge for theorists. 
Therefore, oversimplified models which are still able to 
capture certain qualitative features of the original prob- 
lem have been proposed. The so-called Hubbard model 
_for electrons in a metal and its bosonic counterpart 
assume that the particles only populate a single en- 
ergy band of the periodic potential and that the effective 
interaction has a short range (if not zero range) charac- 
ter. Recent advances in the physics of ultracold atoms 
in optical lattices Q have opened a fascinating frame- 
work in which it is possible to simulate with unprecedent 
accuracy some of the models traditionally used in con- 
densed matter physics. In particular, the transition from 
a superfluid to a Mott insulator of bosons loaded in an 
optical lattice has been successfully observed Q. In an- 
other spectacular experiment [7|] , repulsively bound pairs 
of atoms have been produced, and their main properties 
have been measured. This experiment has stimulated 
renewed interest in the study of few-body effects in dis- 
crete lattices fsl-fl^ , which had been almost only studied 
by Mattis [I^ about twenty years ago. 

Both two-body problems on a one-dimensional tight- 
binding lattice with on-site as well as nearest- 
neighbor [13] interactions can be solved exactly. Sev- 
eral conclusions about them can be obtained (i) 
Bound states can be calculated via a certain polynomial 
equation of different degree depending on the range of 
the interaction, (ii) The scattering states, both symmet- 
ric (bosonic) or antisymmetric (spin-polarized fermionic) 
are well described, asymptotically, by a single phase shift 
which depends again on the range of the potential, (iii) 
The "low-energy" properties of those systems, character- 
ized by the scattering lengths, can always be calculated 
and appear to be rather simple expressions of the respec- 
tive interaction potentials. 



Therefore, the next relevant question to ask is if there 
is a general pattern followed by the solutions and prop- 
erties of the two-body problem when the interactions are 
of arbitrary but finite range. In this article we deal with 
this question and find that, indeed, there is a well de- 
fined pattern for the bound states, and that all scatter- 
ing properties can be calculated very efficiently. For the 
particular, yet very important case of "low-energy" scat- 
tering we show how to obtain the scattering lengths very 
accurately even without knowing the phase shift in gen- 
eral. We treat both identical and distinguishable parti- 
cles which can, of course, have different tunneling rates, 
thanks to a generalization of the center of mass separa- 
tion ansatz for the two-body problem. To illustrate our 
results, we apply them to a model dipolar potential with 
a cutoff at a certain, long enough range. 



II. GENERAL SEPARATION OF THE 
TWO-BODY PROBLEM 

We consider two particles, labeled A and B, in gen- 
eral having different tunneling rates [16] Ja and Jb, and 
interacting via a symmetric two-body potential V{z) = 
V{—z). The reduction of the two-body to a one-body 
problem was first carried out in [l3, [13 • 

The two-body one dimensional discrete Schrodinger 
operator H for the two-body system, acting on a two- 
body wave function u in ^^(Z)(g)^^(Z) is, in first-quantized 
form, 

{Hu){nA,nB) = 

- Ja [u{nA + 1, tib) + u{nA - 1, ns)] 

- Jb [uiuA-, ns + 1) + u{nA, ub - 1)] 

+ V{\nA - nB\)u{nA,nB), (1) 

where ua and n b are the (integer) lattice positions of par- 
ticle A and B, respectively. For simplicity, we set the lat- 
tice spacing s = 1, so that distances, lengths and quasi- 
momenta are dimensionless. For the moment, we do not 
allow V to become infinitely large, that is, < oo 

for all n e Z; this condition can be relaxed by allowing 
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|V(0)| — J> oo and then applying the Bose-Fermi mapping 
theorem (BFMT) ^8j] once the problem is reduced to a 
one-body equation. Moreover, we assume that V is an 
arbitrary finite range potential of range p £ Z, that is, 
V{\n\ > p) = with at least V{p) ^ 0. 
In order to solve this problem exactly we need to trans- 
form the Hamiltonian H to a single particle operator. 
For this purpose, consider the ansatz 

u(n^,nB)=UKWe-'^-^+'^^, (2) 

where R = {nA+nB)/^ and z = UA—riB are, respectively, 
the center of mass and relative coordinates; K is the total 
quasi-momentum and 

tan/3K = 4^-4^ tan (i^/2). (3) 

Note that when Ja — Jr = J, the ansatz ([51 reduces 
to the well known case of identical particles [3-l3- By 
inserting the choice ^ in the Schrodinger equation Hu — 
Eu we arrive at the desired single-particle Hamiltonian 
for each value of the total quasi-momentum K 

{Huk){z) = -\J^''^\[uk{z+1)+uk{z-1)]+V{\z\)uk{z), 

(4) 

where the so-called collective tunneling rate 14] has the 
form 

l^'^^l = ^ J a + JI+'^JaJb cos K. (5) 

Note that the reduced Hamiltonian of Eq. @ is equiv- 
alent to the findings in 113 ■ At this point it is con- 
venient to introduce an adimensional Hamiltonian by di- 
viding it by the collective tunneling, which is equivalent 
to setting J^^^ = 1 (energies become dimensionless) in 
Eq. 0, and rename uk = u for simplicity. We will 
assume this in the subsequent discussions. 

III. BOUND STATES 

We pursue the exact solution for the bound states of 
any two-body system on the lattice with finite range in- 
teractions. Before doing so, we need to define what is 
actually meant by bound state, mathematically, for the 
convenience of the reader. 

We define a bound state of iJ, Eq. ([4]), as any square- 
summable solution u{z) of the discrete time-independent 
Schrodinger equation Hu = Eu with its associated eigen- 
value E lying outside the essential spectrum of H, Ccss = 
[—2,2]. Recall that "outside the essential spectrum" can 
actually mean above 0-S 13 ^'^"^ '^ot only below the 
continuum. 

It is already known that for any finite range potential 
V there exists at least one symmetric bound state [20| : 
it is also known that the maximum number of symmetric 
(antisymmetric) bound states of if is p+1 (p) [21|. Now 
we show rigorously how to calculate all these bound states 



exactly. The formulation of this result is as follows: 
Theorem. 

Let H be the Hamiltonian ^ with V a range-p (< oo) 
potential. Then all bound states u{z) of H have the decay 
property u{z) oc a'^'"'' for \z\ > p, < \a\ < 1; the 
energies of the bound states are given by E = —a — 1/a. 
If u{z) is symmetric then a is a root of a polynomial of 
degree 2p + 1 if p > 1 and, if p = 0, its degree is 2; if 
u{z) is antisymmetric and p > then a is the root of a 
polynomial of degree 2p — 1. 

Proof. Applying the exponential ansatz for u{z) with 
l-z] > p yields immediately 

E = -a- l/a = f{a). (6) 

Since /((-1, 0) U (0, 1)) = (-oo, -2) U (2, oo) and / is 
injective in (—1, 0) U (0, 1), we have that the exponential 
ansatz is the only possible form for the bound states out- 
side the range of V. 

To see that a is a root of a polynomial one shows by in- 
duction, for p > 2, that if w(z) is exponentially decaying, 

then a^u{p — n) = Q^il-ii^) a'^~^u{p — n — 1) = 

*92ri-3^('^)' where Q^™' are polynomials of degree k. For 
symmetric solutions the polynomial equation is then ob- 
tained by setting u(l) = and, for antisymmetric 
solutions, by setting u{0) ~ 0, which proves our state- 
ment. For p — and p = I the result can be proved by 
explicitly obtaining the polynomial equation [9|, [lo| . 

The theorem presented here implies that for any finite 
range potential one has to solve a polynomial equation 
whose degree grows slowly with increasing p. The way of 
obtaining such polynomials is, as can be observed from 
the proof, inductive: we start by setting u{p) = 1 and 
proceed to calculate w(±l) and u(0) by recurrence and 
solve the respective symmetry constrains u(l) — u{—l) 
or u{0) = 0. Certainly if p gets too large it becomes in- 
convenient to get such polynomials for a general potential 
V, and in this case we should obtain the coefficients of 
the polynomial for the given particular potential. 

Consider now the specific choice of the potential 

r ifo< ]z] < 10 

V{z) = I if I2] > 10 (7) 

[ -9.7313 if Izj = 0, 

which corresponds to a dipole-dipole interaction with a 
cutoff at a finite but long range, and where the divergence 
of the potential at 2: = has been substituted by a finite 
value. Note that such dipolar interactions, with a tun- 
able on-site interaction V^(0), can be realized with dipolar 
atoms or molecules in optical lattices [22| - [2^ . We have 
calculated both polynomials P{a) for the symmetric and 
antisymmetric bound states numerically, with their roots 
characterizing the bound states. The results are shown in 
Fig. [TJ For symmetric bound states the polynomial has 
only one root in (—1,1), and therefore only one bound 
state [2^. The polynomial has a root at a = 1, which 
means that it has low-energy resonance. We will discuss 
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FIG. 1: (Color online). Polynomials for symmetric (blue solid 
line) and antisymmetric (red dashed-dotted line) bound states 
with the interaction ((7|, whose roots a characterize the bound 
states with energy E = —a — 1/a. The horizontal line corre- 
sponds to zero ordinate. 



these resonances in Section [V] The polynomial for anti- 
symmetric bound states has also one and only one root 
in (—1,1), in agreement with the discrete Bargmann's 
bound 1261. 



IV. SCATTERING STATES 

After having introduced the first main result of this 
paper, which deals exclusively with bound states, it is 
natural to ask about the exact scattering properties of the 
system. For finite range potentials, the scattering states 
of the Hamiltonian H are asymptotically plane waves, 
that is, for Izj > p we have 



us{z) oc cos(fc|z| -I- 5s), 
ua{z) oc sgn(z) cos(fc|z| -I- 5a), 



(8) 
(9) 



where S and A denote, respectively, symmetric and anti- 
symmetric solutions. Their associated eigenenergies are 
given by the well known tight-binding energy dispersion 
relation [l| 



E = -2cos(fc). 



(10) 



However, a general result concerning the phase shifts 5s 
and 5a does not seem feasible, and it is quite cumbersome 
to obtain them in closed form for long enough ranges. 
One can, however, calculate the phase shifts (and from 
them the exact solution at all z) numerically by recur- 
rence. To this end, we set us{p + 1) = cos(fc(p + 1) + 5s) 
and us{p) = cos(fcp -I- 5s) and analogously for anti- 
symmetric solutions, from equations ([8]) and ([9]). Then 
we calculate us(— 1) and 1*5(1) for symmetric solutions 
with the help of the Schrodinger equation (U) and solve 
ws(— 1) = us(l)- In the case of antisymmetric solutions, 
the relevant equation is wa(0) = 0. We have done so for 
the example potential of Eq. ([7]), as is plotted in Fig. 
[5J There, we clearly observe that the main differences 
between both phase shifts occur at low quasi-momenta 



where the symmetric solution is resonant (see Fig. [Ij; 
looking at slightly higher quasi-momenta already shows 
good agreement between both phase shifts. This means 
that far from fc = fermionization appears rapidly: the 
large on-site interaction V^(0) in (O acts as a hard-core 
at high energies, for which the resonance plays no role (it 
is located at the bottom of the continuum), and therefore 
the symmetric phase shifts are close to the antisymmet- 
ric ( "fermionic" ) . At low quasi-momenta, the resonance 
obviously dominates the asymptotic behavior of the sym- 
metric scattering states. 

In the insets of Fig. [51 we plot the comparison of the 
phase shifts for the potential in Eq. ^ and a model 
range- 1 potential, whose analytic solution is known [lo| . 
The model potential W is chosen so as to be consistent 
with Bargmann's bound [26], and to be resonant for the 
lowest-energy symmetric solution. We obtain [loj 
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W{1) = -^z-^ = -1-19753 

z=l 

W{Q) = -12.125. 



(11) 



The qualitative agreement between the results using V 
or W for the symmetric eigenstates is manifest in Fig. 
[H and, as expected, the differences are most noticeable 
in the high quasi- momentum regime. For antisymmetric 
eigenstates the agreement is very good, even quantita- 
tively, until |fc| ~ 7r/2. The simplified potential (|lip can 
thus be used as a good approximation for the interaction 
([7| in the problem of many (spin-polarized) fermions or 
hard-core bosons on a one dimensional lattice at low en- 
ergies (around the ground state) and low filling (typically 
much smaller than half the number of lattice sites) . Such 
a problem can then be solved exactly by means of the 
Bethe ansatz 12711. 



V. SCATTERING LENGTHS AND 
ZERO-ENERGY RESONANCES 

A. Low-energy scattering 

The "low-energy" (fc — ^ 0, tt) scattering properties of 
the two-body system can be understood via a simple, yet 
exact, calculation of the scattering lengths. Indeed, the 
solution of the time-independent Schrodinger equation 
when fc (fc tt) has an energy E = —2 {E ^ +2), 
and has the asymptotic {\z\ > p) behavior 



usiz) = (Tl)' 



p-ag 



ua{z) = sgn(z)(Tl)' 



(12) 
(13) 



where a~ {a{) is the scattering length at fc -> (fc — tt), 
i = S,A. It must be noted that, in the case of the 
lattice, there are four different scattering lengths, two 
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FIG. 2; (Color online). The calculated phase shifts (tan(5) 
with 5 — 5s or (5 = (5a) for symmetric (blue circles) and anti- 
symmetric (red squares) scattering wave functions, Eqs. ((S)) 
and (|9]), as functions of the relative quasi- momentum k, for 
the potential 0. Left inset: comparison of the symmetric 
phase shift (blue circles) with the one obtained with a model 
range- 1 potential (red triangles), Eq. Hill) . Right inset: an- 
tisymmetric phase shift (red squares) compared to the result 
with the model range- 1 potential (blue triangles). The axes of 
the insets have the same meaning as those of the main figure. 



for "bosons" (symmetric solutions) and two for spin- 
polarized "fermions" (antisymmetric solutions). In order 
to calculate the scattering lengths we proceed as follows 
: using the recurrence relation from z = p by setting 
u,{p + 1) = (Tiy+H^ + 1/(P - a^)], u,{p) = {Tiy 
and E = E± — ±2, the scattering lengths for the 
symmetric states are obtained by solving the equation 
{V{0) - E±)us{0) - 2us{l) = (see proof of the theo- 
rem), while for the antisymmetric states the equation to 
solve is ua{0) = 0. It is remarkable that the resulting 
equations for the scattering lengths as functions of the 
potential can be cast as linear in a^, that is, are of the 
form soa^ + 6o = with sq and bo real constants which 
depend on V{z). In fact, this is an alternative way of 
defining the four lattice scattering lengths, totally equiv- 
alent to the definition a"^ ~ — Um/j_i.^_o dkS [1^, with the 
advantage of not needing to know the phase shift explic- 
itly. It must be noted at this point that, strictly speaking, 
scattering lengths are unique of one-dimensional lattices 
since the radial symmetry is lost in dimensions d > 1. 

As an example, we have calculated for the dipolar 
potential with a cutoff as V{0 < \z\ < p) — —l/\z\'^, 
V{\z\ > p) = 0, and leaving V{Q) as a free parameter. 
The results arc shown in Fig. [3] for a range p = 10, 
where there is a resonance clearly marked at the point 
were the scattering length diverges. 

The divergence of one of the scattering lengths can 
happen for different values of the total quasi-momentum 
K [23] ■ In the simplest case of a zero range interaction 
with V(0) = U, the system is known to have no "zero- 
energy" resonances 0]. For longer ranges, already start- 
ing with p = 1 these resonances can occur. With 
the method outlined in this work we are able to predict 
when, for a given range-p potential with one or more 
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FIG. 3: (Color online). Scattering length (see text) as a 
function of the free parameter V^(O), with V{0 < \z\ < 10) = 
-l/\zf and V{\z\ > 10) = 0. 



free parameters {V{zi), V{z2), . . . , V{zn)}, there is such 
a resonance. To do so, one sets u{z > p) = (t1)^ and it- 
erates recursively as has been explained, and then solves 
the resulting equation for symmetric and antisymmetric 
states, getting a relation among the free parameters of 
the potential to have a resonance. 

We consider again the example of Fig. |3l As we have al- 
ready noted, the system can admit one resonance at the 
bottom of the continuum for the symmetric states. From 
Fig. 131 the approximate location of the resonance can be 
inferred. However, once the scattering length starts to di- 
verge, an accurate location of the resonance grafically is 
very hard (if not impossible), especially if the resonance 
is very sharp as a function of V{0). With our method, 
we are able to locate the resonance very precisely, ob- 
taining for the example we are dealing with a value of 
V{0) — —9.7313 [30|. This is exactly the value chosen in 
the previous sections to match this resonance. 



B. Near- resonant bound states 

It is well known that the binding energy and size of a 
near-resonant bound state (NRBS) are closely related to 
the (large) scattering length ^ 0. On a one dimen- 
sional lattice, the relation between the binding energy E'f, 
and the scattering length in the effective mass approxi- 
mation is given by 



\Eh 



2|m±|(a 



(14) 



where 



=F/i^/2 is the effective mass of the pair with 



an obvious notation. Note that Eq. (1141) is also valid in 
the context of Feshbach resonances |23| . 

When the scattering length is large, the parameter a 
characterizing a NRBS is close to ±1 (for E = —2 and 
E — 2, respectively). By writing \a\ — 1 — |^|, it is 
not difficult to show that the size of the bound state, 
for p > 1, behaves as i* = (|z|) - 1/(2|^|) for |^| < 1. 
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Mathematically, one has 



lim \m* ^ lim 1^ 



{u\ \z\ \u) 



(15) 



Since for small |^|, \Ei,\ we have that the relation 

between the size of the NRBS and its associated scatter- 
ing length is given by 



e* 



2^2 



2 ' 



(16) 



showing that the size of the NRBS grows linearly with 
the scattering length, as we expected. 



C. The number of bound states 

The scattering lengths are very useful quantities in the 
sense that knowing their precise values implies knowing 
the total number of bound states with energies lying be- 
low or above the continuum. For this purpose, we use 
the discrete analog of Sturm oscillation theory [2l|. In 
simple terms, oscillation theory states that the number of 
nodes [si]] of the zero-energy (E = = — 2) symmetric 
(antisymmetric) solution u of Hu — Eu in Z_|_ is exactly 
the number of symmetric (antisymmetric) eigenstates be- 
low E- = —2. Since any state with energy below E- is 
a bound state, the number of nodes of u is the number 
of bound states below the continuum. To see how many 
bound states there are with energies above the contin- 
uum, we make use of the transformation G of Appendix 
1X1 Eq. (|A2[) . or, equivalently, count the number of miss- 
ing nodes. 

We apply oscillation theory now to our example with 
the potential of Eq. (O leaving again V{Q) as a free pa- 
rameter. We calculate (not shown) the symmetric zero- 
energy solution for a given scattering length with 
the methods introduced in this section and see that for 
V{Q) < —9.7313 there are two bound states with energies 
below E- and for V{0) > —9.7313 there is exactly one 
bound state below the continuum. 



VI. A GENERALIZATION 

More general Hamiltonians with exchange operators 
appear when dealing with the problem of one "free" bo- 
son and a bound pair [11]. In such case, the effective 
particles are distinguishable, there is a hardcore on-site 
interaction, an effective range- 1 potential and a first or- 
der exchange term. We assume now that we have the 
following general one-body Schrodinger operator 

{H,^u){z) - - [u{z + 1) + u{z - 1)] + V{\z\)uiz) 

+ n{\z\){Pu){z), (17) 

where V^(O) can be finite or infinite and where P is the 
discrete parity operator. We further assume that ri(|2:|) 



has a finite range pcx with no on-site exchange, J7(0) = 0, 
since it can be included in V{0). 

Obviously [P, Hex] ~ 0, and therefore we can look for 
symmetric and antisymmetric solutions. However, the 
Hamiltonian does not commute with the exchange oper- 
ator QP. With the parity as a good quantum number it 
is straightforward to generalize the theorem of Section HIIl 
to include exchange. To see this, take pm = max(p, pcx)- 
If u(z) is symmetric, the exchange shifts the potential to 
F(|z|) -I- ri(|2:|), while if u{z) is antisymmetric it shifts 
the potential to l^d^l) — f^d^:]). Therefore, obtaining 
the bound states of H^x reduces again to a polynomial 
equation of degree 2pM ± 1, and all the results of our 
theorem apply by changing p by pM and V by V ±rt. 
However, it is no longer true that a hardcore condition 
|y(0)| — ?> oo maps "bosons" onto "fermions" (symmetric 
onto antisymmetric solutions). Indeed, the non-trivial 
dependence of Hex on the parity of the eigenstates makes 
it possible to have states above as well as below its con- 
tinuum even if V and f2 have both the same definite sign. 
This is the fact that makes Hex violate the hypotheses of 
the BFMT |18|], and it explains the appearance of exotic 
three-body bound states on a ID lattice [TTj . 



VII. CONCLUSIONS 

In this paper, we have shown how the exact wave func- 
tions and energies of any bound state of two particles on 
a one-dimensional tight-binding lattice can be calculated 
by solving a polynomial equation whose order increases 
slowly with increasing range of the two-body interaction 
potential, which can also include parity-dependent terms, 
such as effective particle exchange. We have shown that 
the calculation of the exact scattering states is possible, 
and simple. We have also shown how the zero-energy 
resonances associated with the entry or exit of a bound 
state can be trivially and exactly located, and related the 
scattering lengths to the number of bound states above or 
below the continuum by making reference to the discrete 
version of Sturm theory. 

There are, on the other hand, many open problems 
in the physics of few particles on a lattice. At the two- 
body level, there is still no general prescription for the 
calculation of bound and scattering states on two- and 
three-dimensional lattices with arbitrary finite range in- 
teractions. For the three-body case, the Efimov effect 
|32l] , which was shown to appear on a three-dimensional 
(3D) simple cubic lattice [l^, has not been quantitatively 
examined yet; the lattice Efimov effect should depend 
largely on the total quasi-momcntum of the system and, 
moreover, there should appear new kinds of exotic three- 
body bound states in 3D with no analog in continuous 
space, as has been shown to be the case in ID 
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would also be very interesting to explore few-body effects 
in other lattice geometries. 
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Appendix A: Equivalence between repulsive and 
attractive potentials 

We state a general result for an A''-body system on a 
hypercubic lattice in any dimension which, although bc;- 
ing usually implicitly assumed, is useful and important 
to keep in mind, specially when dealing with purely at- 
tractive or repulsive two-body interactions. 

Let H be the following second-quantized Hamiltonian 

H=-jY.ala^ + F{{N}), (Al) 

(m,n> 

where J is the single-particle tunneling rate, (m, n) de- 



notes that the sum runs only through nearest neighbors, 
oj, (fin) is the creation (annihilation) operator of a single 
particle (boson or fermion) at site n and where F is an 
arbitrary analytic function of the number operators at 
each site n, iVn = d\^an. 

If n = (ni, 712, • • • , nd) is a point of a rf-dimensional hy- 
percubic lattice, wc define the unitary operation G so 
that G = & = , with the actions 

Ga], = (-l)^»--a],. (A2) 

We easily see that H — 2F = —GHG~^ is unitarily equiv- 
alent to —H. This implies that the spectrum of a Hamil- 
tonian containing the potentials included in F is obtained 
by changing the sign of every point in the spectrum of the 
corresponding Hamiltonian having F replaced by —F. 
In the case of F containing purely repulsive (or attrac- 
tive) two-body interactions, this result implies the formal 
equivalence between attractive and repulsive potentials. 



[1] N.J. Ashcroft and and N.D. Mermin, Solid State Physics 

(International Thomson Publishing, New York, 1976). 
[2] J. Hubbard, Proc. Roy. Soc. A 276 238 (1963). 
[3] M.P.A. Fisher, P.B. Weichman, G. Grinstein, and 

D.S. Fisher, Phys. Rev. B 40, 546 (1989). 
[4] D. Jaksch et al., Phys. Rev. Lett. 81, 3108 (1998). 
[5] O. Morsch and M. Oberthaler, Rev. Mod. Phys. 78, 179 

(2006); I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. 

Phys. 80, 885 (2008); M. Lewenstein et al, Adv. Phys. 

56, 243 (2007). 
[6] M. Grcincr et. al., Nature 415, 39 (2002). 
[7] K. Winkler et al.. Nature 441, 853 (2006). 
[8] R. Piil and K. M0lmer, Phys. Rev. A 76, 023607 (2007). 
[9] M. Valiente and D. Petrosyan, J. Phys. B 41, 161002 

(2008) . 

[10] M. Valiente and D. Petrosyan, J. Phys. B 42, 121001 

(2009) . 

[11] M. Valiente, D. Petrosyan and A. Saenz, Phys. Rev. A 

81, 011601(R) (2010). 
[12] C. Weiss and H.-P. Breuer, Phys. Rev. A 79, 023608 

(2009). 

[13] L. Jin, B. Chen and Z. Song, Phys. Rev. A 79, 032108 

(2009). 

[14] R. T. Piil, N. Nygaard, and K. M0lmer, Phys. Rev. A 

78, 033611 (2008). 
[15] D. C. Mattis, Rev. Mod. Phys. 58, 2, 361 (1986). 
[16] In continuous space this is equivalent to particles with 

different masses. 
[17] J.-P. Martikainen, Phys. Rev. A 78, 035602 (2008). 
[18] M.D. Girardeau, J. Math. Phys. 1, 516 (1960). 
[19] D. Petrosyan et al., Phys. Rev. A 76, 033606 (2007). 



[20] D. Damanik, D. Hundertmark, R. Killip and B. Simon, 
Commun. Math. Phys. 238, 545 (2003). 

[21] G. Teschl, Jacobi operators and completely integrable 
nonlinear lattices (Mathematical surveys and mono- 
graphs, American Mathematical Society, 2000). 

[22] C. Menotti, C. Trefzger and M. Lewenstein, Phys. Rev. 
Lett. 98, 235301 (2007). 

[23] C. Trefzger, C. Menotti and M. Lewenstein, Phys. Rev. 
A 78, 043604 (2008). 

[24] Th. Lahayc, C. Menotti, L. Santos, M. Lewenstein and 
T. Pfau, Rep. Prog. Phys. 72, 126401 (2009). 

[25] Note that since V{z) < for all z € Z, there can be no 
roots for —1 < a < 0. 

[26] D. Hundertmark and B. Simon, J. Approx. Theory 118, 
106 (2002). 

[27] M. Karbach and G. Miiller, Comp. in Phys. 11, 36 (1997). 
[28] N. Nygaard, R. Piil and K. M0lmer, Phys. Rev. A 78, 
023617 (2008). 

[29] Note that we have normalized the Hamiltonian as H = 

H/|jW]. 

[30] If the calculation is implemented with a symbolic pack- 
age, the value of V'(O) at which the resonance occurs can 
be calculated exactly (with no machine precission limit). 
In our example, V{0) is a rational number whose (large) 
numerator and denominator can be calculated exactly 
this way. 

[31] On the lattice, a function / is said to have a node between 

n and n -I- 1 iff /(n)/(n -h 1) < 0. 
[32] V.N. Eflmov, Phys. Lett. B 33, 563 (1970). 



